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Abstract 

The angular and frequency correlation functions of the transmission coeffi- 
cient for light propagation through a strongly scattering amplifying medium 
are considered. It is found that just as in the case of an elastic scattering 
medium the correlation function consists of three terms. However, the struc- 
ture of the terms is rather different. Angular correlation has a power-law 
decay and exhibits oscillations. There is no "memory effect" as in the case of 
an elastic medium. Interaction between diffusion modes is strongly enhanced 
near the lasing threshold. Frequency correlation scale decreases close to the 
lasing threshold. 

We also consider time correlations of the transmission in the case of non- 
stationary inhomogeneities. We find short- and long-range time correlations. 
The scale of the short-range correlation decreases, while the long-range cor- 
relation scale becomes infinite near the threshold. 
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I. INTRODUCTION 



In recent years novel findings have been uncovered in speckle patterns of waves coher- 
ently transmitted through a disordered medium. The most interesting among them are 
intensity fluctuations with long-range correlation including power-law and "infinite-range" 
correlations. The latter are related to the "universal conductance fluctuations" in disordered 
metals. The theory of speckle arising due to the multiple scattering was developed in |T] [TJ . 
Experimental testing was done for microwave || as well as for visible light H]. 

Recently the study of light propagation in disordered amplifying media has begun. This 
problem is different from a conventional radiation diffusion in a disordered medium because 
here scattering also provides a feedback for light generation and, therefore,upon a certain 
condition an amplifying disordered medium becomes a generator of coherent radiation [0. 
In H authors observed laser action in a dye dispersed in a weakly scattering medium. 
Amplification allows for probing the longest paths that light travels in the system. This 
brings a modification of coherent backscattering || which was observed in [|Kj] . It would be 
natural to carry on these studies by an observation of speckle correlations of light scattered 
from an amplifying medium. 

In this paper we consider correlations in intensity of light transmitted through a dis- 
ordered amplifying medium. Some statistical properties of such a medium were discussed 
TT|Jl^| . Statistical properties of lasing states [|ll|] and transmission distribution func- 



m 



tion [y| were considered. 

Here we calculate the frequency and angular correlation functions of the transmission co- 
efficient of a diffusive amplifying medium. We approach the problem in a way that resembles 
that of H0. 

We show that near the lasing threshold the correlation function consists of three terms 
just as in the case of an elastic medium. However, the structure of the terms is different. 
The first term ( C^ 1 ' ) which describes the contribution of noninteracting diffusion modes 
has a power-law decay . There is no " memory effect" : the term depends separately on the 
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momentum differences of the incoming and outgoing waves. 

The contributions from interaction of diffusion modes are more singular near the lasing 
threshold compared to The term has the same dependence on the incoming and 

outgoing momentum differences as the The "infinite-range" correlation observed in an 
elastic medium is also present in our case, but it has a much larger relative magnitude. 

In frequency dependence of the correlation function near the threshold there appears a 
correlation scale which decreases close to the lasing threshold. The dependence is the same 
in all three terms. 

We also calculated time correlation function in the case of a nonstationary dielectric 
function fluctuations, caused by freely moving point scatterers. We obtain short- and long- 
range time correlations. 

Our paper is organized as follows. In section 2 we introduce the basic quantities which 
we use in our calculations such as the average Green's function of the wave equation and 
the field-field correlation function. In section 3 we define the transmission coefficient and 
calculate the average value of it. We also discuss here the weak localization correction to 
the diffusion coefficient. In section 4 we calculate the angular and frequency correlation 
functions of the transmission coefficient and time correlation in the case of nonstationary 
inhomogeneities. 



We consider the propagation of electromagnetic waves in a weakly disordered amplifying 
medium. We start with a scalar linear wave equation derived from Maxwell's equations: 



Here k = ^ is the wavenumber, e(r) — the complex dielectric function. It is of the form 
e = e' + ie", where the average part: 



II. FIELD-FIELD CORRELATION FUNCTION. 



[V 2 + k 2 e(r)]E = 



(1) 




(2) 
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= 2mn l - Nv y 7 

mc^o - ^o) + 7 

describes the response of oscillating dipoles of frequency w and mass m. N v and N L are 
the densities of the dipoles in the excited and low states correspondingly ( Njj > Nl ). 7" 1 
corresponds to the lifetime of the excited state. 

We will consider the propagation of electromagnetic waves with frequencies oj in an 
interval near the atomic frequency much smaller than broadening 7 so that we can neglect 
the dispersion of the dielectric function. Thus we consider e" as a negative constant and 
suppose that e' « 1. 

Now let some defects or impurities cause fluctuations of the real part of the dielectric 
function. Then: 

e'(f) = 1 + 5e(r) (4) 
where <5e(r) is a random function with the following properties: 

(Se(r)) = 

and 

k\5e{r)5e{f')) = A5(r - r') (5) 

Here angular brackets denote averaging over realizations of the disorder. 
A is simply related to the mean free path of the radiation l : 

4:71 

A = f (6) 

In the lowest order of the disorder parameter (/c/o)" 1 <C 1 the average Green's function 
of the equation (1) is: 

D(r) = — - — exp (ikr -\ (7) 

where 7 = 7- — fc|e"| = 7- — 7^—. 

t tQ t0 tamp 
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The length l amp = ^ott is the amplification length. 

The properties of light energy transport in a disordered medium are determined by the 
field-field correlation function: 

r^fi,^) = (E^n) E:,(f 2 )) (8) 

Here E u is the field of frequency uj. 

Fu;uj(r,r) equals (Iu(r)) — the average intensity at a point r of the field of frequency uj. 

In the leading order of the disorder parameter the field-field correlation function is given 
by the sum of the ladder diagrams: 



ro*/(ri,f 2 ) = — / dr'dr"D(f x - r')D\r 2 - f')P^(r',r")E n u(r")KAr") (9) 
to J 



Here E w {r) is the nonscattered field at a point r and P wu> ' satisfies the equation: 

(dV 2 + — + i 8 J] P^(r) =--5(f) (10) 

where 8uj = uj — uj' , r amp = l am p/c, r = l/c and D = lc/3 is a diffusion coefficient. 
Let us also introduce the light energy current density: 

-* %c — * — * 

U?) = ^(^k - E ^ E -) ( n ) 

It can be easily shown using (9) (see J3j) that: 

(f u (f)) = -DV(L(r,t)} (12) 
Taking into account a conservation law: 

+ ? {m = M2> (13) 

'Tamp 

one obtains a diffusion equation for the average intensity: 

DV\W,t)) + = ^« (14) 

Tamp (ft 

Equations (12) — (14) describe the propagation of radiation in a medium with dimensions 
which are larger than the mean free path Iq. We consider a sample in the form of a slab 



with dimensions L x , L y and L z , where the sides z = 0,L Z are open and the other side 
surfaces are perfectly reflecting. We suppose that L x ~ L y ~ L 2 . The slab has a critical 



width L cr = Ti J DT amp . If L z > L cr the system becomes a multimode generator 0,0]. We 
will measure the size of the system in the units of the critical size: 



Mi - A ) ( 15 ) 



We suppose that A <C 1. 



The case of L z ^> Iq corresponds to the experimental setup of ||10|| . The opposite limit 
L z <C lo, L x , L y was realized in ||. Critical conditions for this case are very different from (15) 
and the critical value of amplification length is of the order of the mean free path |13[ . The 
total phase diagram amplifier-generator is shown on fig. 1. More discussion of the critical 
conditions can be found in fL4| . Calculations in the case L z < l are beyond the scope of 



our paper. 

III. THE AVERAGE TRANSMISSION COEFFICIENT 

Consider the transmission of monochromatic light through the slab described in the 
previous section. 

We define the transmission coefficient of the light of frequency u incident at an angle a 
on the plane z = and transmitted at an angle b as the ratio of the emergent flux per unit 
solid angle around the direction % to the incident flux. 

From (9) and (11) one has for the transmission coefficient: 



Ic cos 2 b r .o d , , _.. , — Z 3 cos 2 b f , ,„ , d d 

(16) 



where Jo is the incident intensity. 

The corresponding to the equation (16) diagram is shown on fig. 2. 
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For the equation (10) for function P^'if^r') we assume the following boundary condi- 
tions: P^uj' is zero on the boundaries z = 0, L z and the current through the side surfaces is 
zero. Then near the threshold one can calculate the function P WW '(f, r') as: 

P ,(r r>) - 3Lcr Sm ^ z)sm &) (m 
^ [r ' r) ~lW x L y A(l-zf) [ } 

where Q = 2A/r amp . In (17) we neglect contributions which are nonsingular at A — > 0. 
Hence one obtains: 

3/ 

(T abw ) = j—£ cos 2 b (18) 
Integration over angles gives the average coefficient of the total transmission: 

< T > = < 19 > 

The transmission coefficient contains an additional amplification factor A -1 compared to 
the case of an elastic medium. 

Let us now consider the weak localization correction to the diffusion coefficient in equa- 
tion (14) in the limit A — > 0. Weak localization has been extensively studied in physics of 
disordered conductors ( see for review \TE\ ) and we may use the well known results. Lo- 



calization correction to the diffusion constant is given by maximally crossed diagrams. The 
analytical expression is: 

D(5u) = D M - My dfC^(f,f)j (20) 

Here V = L x L y L z is the volume of the slab and cooperon Cw(r, r ') satisfies in the case of 
stationary inhomogeneities eq.(10) and has therefore form (17). Substituting it in (20) one 
obtains: 
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D{5uj) = D 



1 



i TrAfl-i^),/ 



(21) 



Here g = - L L xl ^ yl ^> 1 is a quantity analogous to the dimensionless conductance of a weakly 
disordered metal. Let us note that there is no correction to l/r amp in eq.(14). For A — > 
the correction in (21) coincides with that for electrons in a small closed metallic particle. 



For small Sou, such as 5u ~ ~ 7^7, the correction is of order unity. In our case it 

means that at a time ~ (du)" 1 diffusion approach breaks down and transport is dominated 
by only a few lasing states. At zero 5uj the condition of smallness of the correction is g A > 1 . 



IV. INTENSITY-INTENSITY CORRELATIONS 

Average quantities do not describe the transmission of light in a random system com- 
pletely. In this section we consider fluctuations of the transmission. We study the following 
correlation function: 

Caba'V = ^T abul 5T al y^) (22) 

where 5T ahuJ = T abuJ - (T abuJ ). 
It might be represented as ||: 

C%£ v = C w + C (2) + C (3) (23) 

Contribution describes noninteracting diffusion modes. Interaction between them brings 
contributions and C^ 3 -*. The diagram corresponding to C^~> is shown on fig.3. The 
calculations leading to an expression for are analogous to the calculations of the average 
transmission coefficient. The only modification needed is to take into account a change in 
the phase of the incident or transmitted wave due to a change in the angle of incidence or 
transmission: 

d d 



^) I 6 cos 2 b cos 2 V 



J d 2 r d 2 r' exp (-i5q b r) exp (i5q a r') — — 



Ar,r z= r z , =0 



(24) 



Here 5q a ^ is the difference between the wave vectors of the incident (a) or transmitted(fe) 
waves corresponding to different measurements. 
From (17) and (24) one obtains: 

= {T ahw ){T albl ^)F{5q a )F{5q b ) 1 — 1 (25) 

1 + (§) 



where 

sin 2 (8q x If) sin 2 (5q y ^) 
*y d( l) = —, TT2 7 TT2— 

C^ 2 ^ term can be represented as C^> = C^ 2 \5q a ) + C^ 2 \5qb). Due to reciprocity relations 
both terms are equal for 5q a = 8q\. 

Let us note the following. The scatterers in the disordered medium can be considered 
as sources of light, which are phase or amplitude correlated. Accordingly we can divide the 
terms constituting the speckle correlation function in two groups. Terms C^ 1 ' and C^{5qb) 
correspond to phase correlated, C^ 2 '(Sq a ) and — to amplitude correlated sources. 

To calculate the amplitude correlated terms it is convenient to use the Langevin ap- 
proach p. According to this approach the fluctuating part SJ^ = J w — (J u ) of the light 
energy current can be expressed as: 

5Uf) = -DVSUf) + j extuJ (r) (26) 

where jextuir) is the Langevin random current. Using conservation law: 

WJ w (r) = (27) 

7~amp 

one can write for the fluctuating part of the intensity: 

(DV 2 + —) 6L(f) = Vj extuj (r) (28) 

With the help of (26) and (28) we can obtain an expression for the fluctuating part of the 
total current through the surface z = L z . To calculate the correlation function of the total 
current we need an expression for the correlation function of the Langevin current. 

Near the lasing threshold the Langevin currents correlation functions are given by the 
following expressions: 

2tt/c 2 _ 1 

"n J 



(JU u (fi)jLA^))i = -^(Wi)) 2 Fm i - - r 2 ) (29) 



(jL»{n)jLA?2)h = h^r + ^ + k 3 ) (so) 



K x = SijSin - f 2 ) J df 3 \P^if u f 3 )\ 2 V(Ur 3 ,a)}V(I^(f 3 ,a')} (31) 
K 2 = |P^(f 1 ,f 2 )| 2 V J (^(i?i,a'))V,(/ w (f 2 ,a)) (32) 

K 3 = |P^(n^2)| 2 V i (/.(f 1 ,a))V J (/^(f 2 ,a / )) (33) 

Correlation functions (29) and (30) contribute to C^ 2 \8q a ) and correspondingly. Di- 
agram for correlation function (29) is shown on fig. 4 and diagrams for the three terms in 
(30)— on fig.'s 5a,b,c @. 

Calculating fluctuations of the total current we obtain: 

= CW(8q a ) + C (2 \5q b ) = -i—^^w) L-,[F(6&) + F(6q b )] (34) 

AngA x + ^ 

27 1 

ci3) = iwg*A* {T ^ ){T *™\ + f*y (35) 

Let us now discuss our main results (25), (34) and (35). Note first, that there is no 
"memory effect" in the term — it depends separately on 5q a and 5cfb. The form-factor 
function F(8q) has an oscillating power-law decay . 

The form-factor is exactly the same as in the case of transmission through an elastic 



scattering medium in the form of a long tube [16|. The "memory effect" is also absent in 
such a medium. 

Integrating over the angles of transmission we obtain the mean square fluctuation 



of the total transmission 11 



{{ST)2) = 4^a5 (T>2 (36) 
One sees that the transmission near the lasing threshold is not a self-averaging quantity: 
fluctuations of total transmission grow faster than the average value. 

The C (3) term just as in the case of an elastic medium represents an "infinite-range" 
correlation in momentum space. However, the relative magnitude of this correlation is much 
larger in our case due to the factor A~ 4 . 
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Close to the threshold all terms in the speckle correlation function have the same fre- 
quency dependence. The scale of the frequency correlation is determined by a parameter 
Q = 2A/r amp , which goes to zero at A — > 0. 

Interaction between diffusion modes according to (34) and (35) is described by a param- 
eter 1/gA 2 , contrary to the case of a nonamplifying medium where the parameter is 1/g. 
The increase of interaction can be interpreted as follows. In the case of an elastic disordered 
system transmission probes the path of typical length of order L\jl. In our case amplifica- 
tion compensates for escape of a wave, therefore transmission probes paths of much longer 
lengths and probability of interference increases. 

We should point out that the value of parameter A in our formulas can not be infinitely 
small. A = corresponds to the threshold in the equation for the average intensity. 
Actually the threshold value of A is a random function of frequency, such that its mean 
square fluctuation is of the order of g^ 1 It means that our consideration is valid for 



gA 2 > 1. Note that this condition is more restrictive than the condition of smallness of the 
weak localization correction gA > 1 in eg. (14). 

So far we assumed that the fluctuations of the dielectric function which scatter light are 
stationary. However it is also interesting to consider the correlation function C a t, a b'(t) = 
(<5T a {,(t)(5T a b/(0)) ( T a b(t) is the transmission coefficient T a b measured at time t) in the case 
of nonstationary fluctuations. Note that we will consider correlations with a fixed incident 
channel (a = a'). 



In the simple case of freely moving point scatterers [17. 18] one obtains: 



C^\t) = (T ab )(T ab> )F(5q b )F 1 (t) (37) 



C (2) (*) = — ^ (Tab) (Tab') [F 2 (t) + F(6q b ) F 3 (t)] (38) 



Here 

1 



Flit) 



2 

i -r o^ri 



2t^ fir 
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DC 



F 2 (t) = n 2 / dh dt 2 ex P {-n (ti + t 2 )} F x (h -t 2 + t) 



o 



F 3 (t) = VI 2 J™ dti dt 2 exp j-fi ( 1 + 2r 2 flr ) (*i + - h) 

and ta is the time for a scatterer to move a wavelength. We do not consider the 
term because in the present case it does not contain additional features compared to 
and C (2) . 

To derive the above results one needs to modify the dielectric function fluctuation's 
correlation function (5) to take into account the motion of the scatterers. In the case of 
freely moving scatterers with a characteristic velocity much less than the velocity of light in 



the medium this correlation function averaged over the angle of scattering becomes 117,18 



k 4 (5e(f,t)6e(f',0)) = Af(t)5(f- r') (39) 

where f(t) = (r 2 /t 2 ) [1 — exp(— t 2 /r 2 )} ps 1 — t 2 /2r 2 and t is the time difference variable. 

Then it is easy to show that function P uu i(t) describing the correlation of fields 
(£ , CJ (t)i?*/(0)) measured with a time interval t satisfies the equation: 

(DV 2 + — + i 6oo - P^{r,t) = - - 5(r) (40) 

Therefore near the threshold it is given by an expression: 

P uu >(r,r ,t) = — — (41) 

The further calculations follow the same lines as those in the case of stationary scatterers. 
According to (37) and (38) there exist short-range time correlations Fi(t) and 

A 3/2 



2t 2 Qt 



with a scale t\\/VLt and long-range time correlation F 2 (t) ~ A 3 / 2 exp (—fit). 
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V. CONCLUSIONS 



In conclusion, we have calculated angular, frequency and time correlation functions for 
the transmission of light through a disordered amplifying medium. The correlation function 
in the case of stationary dielectric function fluctuations consists of three terms like in an 
elastic scattering medium. The angular dependence of the correlation function is exactly 
the same as in an elastic medium with tube geometry. This similarity is due to the absence 
in the emitted radiation of memory about the incident radiation in both cases. Frequency 
correlation has a scale Q = 2A/r amp — > at A — * 0. It is the same for all three terms. 

Interaction between diffusion modes is strongly enhanced near the lasing threshold. The 
enhancement of transmission fluctuations is more pronounced than the enhancement of weak 
localization near the threshold. 

In the case of nonstationary dielectric function fluctuations, speckle correlation function 
consists of terms with short- and long-range time correlations. The short-range correlation 
has a scale t\V^t — » at A — > 0. Interaction between diffusion modes is suppressed by the 
motion of the scatterers. 

Although our calculations were performed in the case of all sample dimensions being 
larger than the mean free path, we believe that our results also apply to the case L z < l . 
Probably only the "interaction parameter" (gA 2 ) -1 and definition of parameter Q will be 
modified. 

Experimentally it would be interesting to realize the crossover between the both cases 
in one system. Probably, a good candidate is an amplifying liquid crystal near the phase 
transition, which unduces strong changes of transparency. Recent experiments Jl9| show 
that it is possible to realize the multiple scattering regime in a liquid crystal. 
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FIGURE CAPTIONS 



1. Figure 1. 

Schematic picture of an amplifier-generator phase diagram for a disordered slab of thick- 
ness L z with a homogeneous amplification length l amp - In the region below the curve the 
slab is a random amplifier. Above the curve it is a random generator. 

2. Figure 2. 

Diagram for the diffusion propagator P. Solid lines correspond to the average Green's 
function D, dashed lines denote scattering. 

3. Figure 3. 

The diagrammatic representation of C^. 

4. Figure 4. 

The diagrammatic representation of the correlation function of the Langevin currents 
which are proportional to the average intensity. 

5. Figures 5a,b,c. 

The diagrammatic representation of the correlation function of the Langevin currents 
which are proportional to the average current. We denote the average current by a shaded 
vertex. 
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